Experimental probing of exchange interactions between localized spins in the dilute 
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The sign, magnitude, and range of the exchange couplings between pairs of Mn ions is determined 
for (Ga,Mn)N and (Ga,Mn)N:Si with x < 3%. The samples have been grown by metalorganic vapor 
phase epitaxy and characterized by secondary-ion mass spectroscopy; high-resolution transmission 
electron microscopy with capabilities allowing for chemical analysis, including the annular dark-field 
mode and electron energy loss spectroscopy; high-resolution and synchrotron x-ray diffraction; syn- 
chrotron extended x-ray absorption fine-structure; synchrotron x-ray absorption near-edge structure; 
infra-red optics and electron spin resonance. The results of high resolution magnetic measurements 
and their quantitative interpretation have allowed to verify a series of ab initio predictions on the 
possibility of ferromagnetism in dilute magnetic insulators and to demonstrate that the interaction 
changes from ferromagnetic to antiferromagnetic when the charge state of the Mn ions is reduced 
from 3+ to 2+. 



I. INTRODUCTION 

The decisive role of holes in ordering the localized spins 
in dilute magnetic semiconductors (DMSs) is not only 
well established? 1 - - — but it represents also the basis of 
the functionalities demonstrated for these systems^ In 
view of the fact that most magnetic insulators are either 
antiferromagnets or ferrimagnets, particularly intriguing 
is the question whether ferromagnetism is at all possi- 
ble in dilute magnetic insulators, where carriers remain 
strongly localized on parent impurities or defects^^ Ac- 
tually, a ferromagnetic coupling between localized spins 
was predicted in a series of ab initio works for the model 
system (Ga,Mn)N^ where, as shown in Fig. [I] a large 
value of the ferromagnetic exchange energy J nn was cal- 
culated for the nearest neighbor (nil) Mn pairs. Due 
to the highly localized character of the orbitals in ques- 
tion, the magnitude of J is expected to decay fast with 
the inter-spin distance. Nevertheless, according to recent 
Monte Carlo simulations, the predicted Curie tempera- 
ture Tq is as high as 35 K and 65 K for the Mn cation 
concentration x = 3% and 6%, respectively^ 

As reviewed elsewhere this clear-cut theoretical 
prediction has not been yet verified experimentally. In- 
stead, a diversity of magnetic properties has been re- 
ported. For instance, no indication of ferromagnetic in- 
teractions was detected up to x = 36% in polycrystallinc 
films prepared by ion-assisted deposition^ whereas Tq 
values ranging from 8 K (Ref. [l2l ) up to over 300 K 
(Ref. H) were found for (Ga,Mn)N grown by molecular 
beam epitaxy (MBE). However, the detection of phase 
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FIG. 1: (Color online) Exchange energy J nn for the nearest 
neighbor (nn) coupling from ab initio computations by vari- 
ous authors, as listed in Ref. 0- To obtain J nn , the calculated 
energy differences for antiferromagnetic and ferromagnetic ar- 
rangements of spins within the nn Mn pair of the determined 
exchange energies have been mapped on the classical Heisen- 
berg Hamiltonian. 



nanoparticlcs. It is increasingly clear that a further 
progress in the understanding of this challenging system 
requires a precise control of both the spatial distribution 
and the charge state of the Mn ions<2 

The samples whose properties are discussed in this 
paper have been grown by matalorganic vapor phase 
epitaxy (MOVPE) according to a procedure described 
in Sec. II. As detailed in Sec. Ill, our films have 
been characterized by: secondary-ion mass spectroscopy 
(SIMS); high- resolution (scanning) transmission elec- 
tron microscopy [HR-(S)TEM] with capabilities allow- 
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ing for chemical analysis, including the energy disper- 
sive x-ray spectroscopy (EDS), high angle annular dark- 
field (HAADF) mode, and electron energy loss spec- 
troscopy (EELS); high- resolution and synchrotron x- 
ray diffraction (SXRD); synchrotron extended x-ray ab- 
sorption fine-structure (EXAFS); synchrotron x-ray ab- 
sorption near-edge structure (XANES); infra-red optics 
and electron spin resonance (ESR). This set of probes 
demonstrates the absence of precipitation, reveals a non- 
uniform Mn distribution in the digitally Mn-doped films, 
and shows that the concentration of Mn 2+ ions reaches 
4 x 10 20 cm -3 in (Ga,Mn)N:Si. From four probe conduc- 
tance measurements, the sheet resistance is of the order 
of 10 GQ at 300 K for the (Ga,Mn)N film with the high- 
est Mn content x = 3.1%, where x is the concentration 
of Mn cations. 

By combining this extensive growth and nanochar- 
acterization program with the results of high-precision 
magnetic measurements discussed in Sec. IV we demon- 
strate that the dominant interactions between neighbor 
Mn pairs are ferromagnetic in (Ga,Mn)N. However, ac- 
cording to the data, the coupling is too short-ranged to 
lead to magnetic ordering above 1.85 K in the studied Mn 
concentration range up to x = 3%. Employing a model 
of magnetic susceptibility suitable for wurtzitc (Ga,Mn)N 
at high temperatures, and outlined in the Appendix, we 
evaluate from our experimental results the magnitude of 
the exchange energy for the nearest neighbor ferromag- 
netic coupling. These findings allow to verify the series 
of ab initio predictions summarized in Fig. Q] on the pos- 
sibility of ferromagnetism in dilute magnetic insulators. 
At the same time, we show that the interactions become 
antifcrromagnetic if the Mn charge state is altered by co- 
doping with Si donors, clarifying in this way the array of 
magnetic properties reported for this system. 



Sample 


Label 


Mn x (av) 


Thickness 


Ga 


Mn 


Si 


# 




(%) 


(nm) 








966 


(Ga,Mn)N 


0.5 


470 


5 


490 





1069 


(Ga,5Mn)N 


1.8 


140 


5 


(5-490 





1071 


(Ga,5Mn)N 


2.6 


135 


5 


5-490 





1080 


(Ga,Mn)N 


1.1 


740 


5 


490 





1106 


(Ga,Mn)N 


1.8 


750 


1 


490 





1130 


(Ga,Mn)N 


1.8 


200 


1 


490 





1134 


(Ga,Mn)N:Si 


1.8 


220 


1 


490 


2 


1142 


(Ga,Mn)N 


3.1 


230 


1 


490 





1152 


(Ga,Mn)N:Si 


3.3 


200 


1 


490 


2 


1159 


(Ga,5Mn)N 


1.5 


140 


1 


5-490 





1160 


(Ga,5Mn)N:Si 


2.4 


135 


1 


5-490 


2 


1161 


(Ga,Mn)N:Si 


2.9 


200 


1 


490 


1 


1268 


(Ga,Mn)N:Si 


2.0 


232 


1 


400 


2 


1269 


(Ga,Mn)N:Si 


1.9 


232 


1 


300 


2 


1273 


(Ga,Mn)N:Si 


0.49 


780 


5 


490 


2 


1274 


(Ga,Mn)N 


0.53 


780 


5 


490 






TABLE I: Samples studied in this work with Mn cation con- 
centrations x and x av for the uniformly [(Ga,Mn)N] and dig- 
itally Mn-doped structures [(Ga,5Mn)N], respectively, as de- 
termined by fitting the magnetic model to the data obtained 
by superconducting quantum interference device (SQUID), 
as described in Sec. IV. The total thickness of the Mn-doped 
layer and precursor flow rates (in seem) for the Ga, Mn, and 
Si containing precursors are also given. The samples #1273 
and 1274 have been grown on double-side epi-ready sapphire 
substrates suitable for optical transmission studies. 

(Sec. IV). An additional series of samples has been grown 
onto double-side epi-ready substrates for optical trans- 
mission studies. The Mn concentration x (x av for the 
digital structures), as determined by the near-saturation 
value of the in-plane magnetization M at 50 kOe and 
1.85 K, reaches over 3% in the samples with the highest 
Mn content. 



II. GROWTH METHOD AND STUDIED 
SAMPLES 

In order to increase the Mn concentration in (Ga,Mn)N 
grown on GaN/c-sapphire by MOVPE at a substrate 
temperature of 850° here the flow-rate of the Ga pre- 
cursor (TMGa) is reduced to 1 standard cubic centimeter 
per minute (seem), maintaining the temperature of the 
Mn precursor source (McCp2Mn) at 22°C and its flow 
rate up to 490 seem. In addition to the uniformly doped 
Mn films, we also grow digitally (6) Mn-doped structures, 
in which the Mn and Ga precursors are supplied alter- 
nately with a period ratio up to 8. Furthermore, a series 
of respectively uniformly and digitally Mn-doped sam- 
ples is co-doped with Si at a S1H4 flow rate of either 
1 or 2 seem. The four types of considered samples are 
denoted by (Ga,Mn)N, (Ga,5Mn)N, (Ga,Mn)N:Si, and 
(Ga,<5Mn)N:Si respectively, where with 5Mn we refer to 
the Mn-digitally doped layers. 

In Table I the list of the studied samples is presented, 
whose magnetic properties are reported in Figs. llll and ll2l 



III. NANOCHARACTERIZATION 
A. Crystallinity 

The degree of crystallinity and possible precipitation of 
secondary crystallographic phases in uniformly and dig- 
itally Mn-doped films have been assessed by HR-TEM, 
HR-STEM, HR-XRD, SXRD, and EXAFS. 

Our TEM studies have been carried out for all studied 
films on cross-sectional samples prepared by mechanical 
polishing followed by Ar + ion milling, under a 4° angle 
at 4 kV for less than 2 h. The ion polishing has been 
performed in a Gatan 691 PIPS system. 

The specimens have been investigated in Linz using a 
JEOL 2011 Fast TEM microscope operating at 200 kV 
and equipped with a Gatan CCD camera. The set-up 
is capable of an ultimate point-to-point resolution of 
0.19 nm, with the possibility to image lattice fringes with 
a 0.14 nm resolution. 

As reported in Fig. for both (Ga,Mn)N (x = 3.1%) 
and (Ga,<5Mn)N (x av = 2.6%) samples low resolution 
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FIG. 2: Low-magnification TEM (left panel) and HR-TEM 
(right panel) for (Ga,Mn)N (x = 3.1%, upper panels) and 
(Ga,<5Mn)N (x av = 2.6%, lower panels) without any evidence 
of crystallographic phase separation. 



TEM (left panels) shows no crystallographic phase sep- 
aration and, in particular, no precipitates' segregation 
near the surface or interface. In fact, the HR-TEM im- 
ages (right panels) clearly reveal the atomic positions in 
the lattice and, on the scale displayed, they show a ho- 
mogenous crystal ordering and no signs of precipitation 
within the Mn-doped layers. 

To further verify the crystallographic homogeneity of 
the grown Mn-doped layers, HR-XRD measurements us- 
ing a Materials Research Diffractometer (MRD) and 
SXRD were performed. The HR-XRD experiments have 
been carried out with a Panalytical X'Pert PRO MRD 
in Linz at the photon energy of the Cu K a i radiation 
(8 kcV) using a hybrid monochromator with a 0.25° slit 
for the incident optics and a pixel detector with an active 
length of 1 mm (19 channels) in the diffracted beam op- 
tics. The SXRD experiments have been performed at the 
beamline BM20 (Rossendorf Beam Line) of the European 
Synchrotron Radiation Facility (ESRF) in Grenoble us- 
ing a point detector and 0.5 mm slits in front of the beam 
at an energy of 10 keV. Radial ui-29 scans of the GaN 
(002) to the GaN (004) diffraction peak do not show any 
trace of secondary phases. In Fig. [3] the radial scans ac- 
quired with both techniques for the same set of samples 
are shown for comparison, with the MRD scans in the 
upper panel and the SXRD in the lower one. The sharp 
(002) and (004) diffraction peaks of GaN and the (006) 
diffraction of the sapphire substrate are observed for all 
samples. The measurements carried out with the high 
energy monochromatic beam at the synchrotron show 
a better signal-to-noise ratio and less diffuse scattering, 
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FIG. 3: (Color online) XRD radial scans of (Ga,Mn)N and 
(Ga,<5Mn)N: (a) HR-XRD MRD, (b) SXRD collected at the 
BM20 of the ESRF. 



when compared to those performed with the MRD, but 
in both cases no trace of crystallographic precipitation 
can be observed. 

The narrow full- width at-half-maxima (FWHM) of the 
GaN (002) and (004) peaks (with values between 240 and 
290 arcsec) indicate the high crystallinity of the layers. 
From the GaN symmetric sharp diffraction peak from 
the Mn-doped samples - comparable to the one from the 
GaN reference - we have hints that Mn-doping does not 
affect critically the dislocation density, as confirmed by 
the HR-TEM analysis on the same samples. 

The XAFS measurements at the Mn K-edge (6539 eV) 
have been carried out at the GILDA Italian collaborat- 
ing research group beam-line (BM08) of the ESRF in 
Grenoble, according to the experimental procedure de- 
tailed previouslyJii£ Both EXAFS and XANES regions 
of the collected spectra have been analyzed. 

Three representative samples have been studied: 
(Ga,<5Mn)N (x av = 2.6%), (Ga,Mn)N (x = 3.1%) and 
(Ga,Mn)N:Si (a; = 3.3%). The collected data, with the 
polarization vector parallel to the c-axis, and the rela- 
tive fits resulting from the EXAFS analysis are shown 
in Fig. The structural model employed consists of 
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FIG. 4: (Color online) Magnitude of the Fourier-transformed 
fc 2 -weighted EXAFS data (inset) in the range fcmin-fenax with 
relative fits (described in the text) in the range i? m in-.Rmax for 
representative samples: (Ga,<5Mn)N (a; a v = 2.6%), (Ga,Mn)N 
(x = 3.1%) and (Ga,Mn)N:Si (x = 3.3%) 



one Mn atom substituting Ga in a GaN wurtzite crystal 
plus a Mn-Mn contribution taken from the MnN crystal 
structure^ with a fitted coordination number (Nmh-Mu) 
in order to account for possible Mn clusters. The quan- 
titative results are reported in Table [TT| and they are 
equivalent within the error bars for all the samples. The 
agreement with experimental data is good up to several 
coordination shells - as seen in Fig. @] - demonstrating 
the substitutional incorporation of Mn. The presence 
of an additional Mn-Mn coordination can be excluded 
within an uncertainty of 0.4 neighbors. Considering that 
the typical metal-metal coordination numbers are around 
10 (12 Mn neighbors in MnN, 8 in Mn 3 N 2 ) we can ex- 
clude the presence of secondary phases with a precision 
of about 5%. In addition, the possible Mn incorporation 
in interstitial sites (tetrahedral and octahedral) has been 
investigated, statistically not improving the fit quality. 
These results are equivalent to the XAFS structural anal- 
ysis on our (Ga,Mn)N samples at lower concentrations 

In conclusion, TEM, HR-(S)TEM, HR-XRD, SXRD, 
and EXAFS studies demonstrate a single-crystal char- 
acter of the considered films with no traces, down to 
the atomic scale, of precipitations (crystallographic phase 
separation). This is in contrast with the case of some 
(Ga,Mn)N films grown by MBE (Refs. EHJ) as well as 
with (Ga,Fe)N layers obtained by MOVPE (Refs. GJIH- 
|24| ) , for which crystallographic phase separations were de- 
tected by employing the same nanocharacterization tools. 



B. Mn distribution 

The Mn distribution along the growth direction and 
the Si distribution in the co-doped films have been 
evaluated in Warsaw via SIMS, calibrated by Mn im- 
planted GaN, providing the absolute concentration of Mn 



TABLE II: Quantitative results of the EXAFS analysis. For 
each sample, the fitted parameters are: the common ampli- 
tude (So), the average bond distances from the central Mn to 
the 4 N nearest neighbors (Rmh-n) and 12 Ga next nearest 
neighbors (7?Mn-Ga) of the wurtzite structure plus a common 
expansion parameter for higher coordination shells (AR) and 
the coordination number (Nmi-Md) f° r the Mn-Mn bond dis- 
tance at 2.98(2) A. The Debye- Waller parameters attest all 
below 8(2) x 1(T 3 A" 2 and a correlated Debye model 2 ^ with 
a temperature of 500(50) K is used for the GaN multiple scat- 
tering contributions. Error bars are reported on the last digit. 



Sample 


^(av) 


(%) 


flMn-N 

(A) 


-RMn-Ga 

(A) 


AR 

(%) 


NMn-Mn 


(Ga,5Mn)N 

(Ga,Mn)N 

(Ga,Mn)N:Si 


2.6 
3.1 
3.3 


0.90(5) 
0.94(5) 
0.94(5) 


1.95(3) 
1.94(3) 
1.96(3) 


3.20(2) 
3.19(2) 
3.19(2) 


0.1(2) 
0.1(2) 
0.1(2) 


0.4(4) 

0.2(4) 
0.0(4) 



atoms with an accuracy of about a factor 2. The SIMS 
depth profiles reported previously^ for the samples with 
x < 1% and the one shown in Fig. [5] for the present 
films indicate that the Mn distribution is uniform over 
the Mn-doped region and that the interface between the 
(Ga,Mn)N ovcrlayer and the GaN buffer layer is sharp. 




Depth (urn) 

FIG. 5: (Color online) SIMS depth profiles of Mn and Si in the 
(Ga,Mn)N:Si film with x — 1.9%. The residual contamination 
by H, O, Fe, and Mg is also shown. 

This is confirmed by chemical analysis performed in 
Linz with an Oxford Inca EDS system, which — with the 
sensitivity of 0.1% at. — does not provide any evidence for 
Mn diffusion into the nominally undoped GaN buffer. 

For further characterization of the Mn doping, 
HAADF-STEM and EELS measurements have been per- 
formed - and reported in Fig. [5]- employing a FEI Tecnai 
F 20 200kV transmission electron microscope in Graz. 
The STEM images and EELS spectra could be recorded 
with an upgraded spectrometer with an adapted STEM 
detection geometry, optimized for Z-contrast and en- 
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hanced spectral collection EELS efficiency! 25 ' 26 In order 
to improve the signal to noise ratio at the high spatial res- 
olution in question as well as to keep exposure times low 
for minimal sample drift, the EELS point spectra have 
been taken as an integral sum over the energy. Accord- 
ingly, a possible fine structure located within the first 
20-30 eV around the ionization threshold energy of an 
edge has been averaged out. 

The HAADF-STEM on (Ga,Mn)N (x = 3.1%) and re- 
ported in Fig. EJa) reveals a clear change in intensity, 
when going from the substrate into the Mn-doped layer 
(from left to right), while within the doped layer no chem- 
ical contrast could be detected. The EELS spectra given 
in Fig. [l^b) evidence the presence of Mn only in the nom- 
inally Mn-doped layer. 

In contrast, remarkably, as shown in Figs. [SJc) 
and life), HAADF-STEM observations on the digital 
(Ga,<5Mn)N (x av = 2.6%) reveal intensity modulations. 
Additional EELS point spectra collected in Fig. EJd) 
show the signal differences of the Mn L23 edge between 
the substrate and the intensity modulated lines. 

In conclusion, the clement specific analysis demon- 
strates a spatially homogeneous Mn distribution over the 
volume of the uniformly Mn-doped (Ga,Mn)N films, with 
no segregation towards the surface, interface or buffer re- 
gions. In contrast, in the case of the digitally Mn-doped 
films, nano-scale density modulation with a period im- 
posed by the growth conditions has been detected, mean- 
ing that in these <5Mn samples the local Mn concentration 
fluctuates between lower and higher x values around the 
x av determined by SQUID magnetometry. As discussed 
in Sec. IV, such a non-random distribution of Mn ions 
increases the apparent Curie constant, particularly if the 
system is close to a ferromagnetic instability. 
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C. Concentrations of Si and Mn 2+ ions 

The incorporation and a uniform distribution of Si im- 
purities in (Ga,Mn)N:Si layers is evidenced by the SIMS 
result displayed in Fig. [5j From the same measurements 
on all the considered samples, the Si concentration is 
found to be of the order of 10 20 cm -3 for the Si and 
Ga precursor flow rates 1 or 2 and 1 seem, respectively. 

Following recent works of the Mn charge state in 
GaN;22 the effect of co-doping by Si donors on the Mn 
charge state has been quantitatively assessed by exam- 
ining the magnitude of the intra-ion optical absorption, 
which occurs at E a « 1.4 eV for Mn 3+ ions in GaN.^^L 
Optical investigations have been performed in Warsaw 
and in Linz for two series of films, abridged in Table [TTTl 
with the Mn concentrations x w 1.8 and 0.5%, respec- 
tively, and different Si content. The series with the lower 
Mn concentration has been designed for transmission and 
ESR measurements with the layers deposited on double- 
side epi-ready sapphire substrates. 

As it can be seen in Figs. [7] and [SJ the Si dop- 
ing quenches the intra-ion absorption, specific to Ga- 



FIG. 6: (Color online) (a) HAADF-STEM image of 
(Ga,Mn)N (x = 3.1%)- change in intensity (chemical con- 
trast), when going from the substrate into the doped layer 
(from left to right), (b) EELS spectra for points (GaN 
buffer layer) and 1 (nominally Mn doped layer): evidence of 
the presence of Mn only in the Mn-doped layer (right side of 
the image), (c-f) Determination of the Mn distribution in the 
digitally Mn-doped sample (Ga,<5Mn)N (a; av = 2.6%, 160 Mn 
periods, layer thickness 135 nm). (c,e) HAADF-STEM scans 
in the Mn-doped region giving modulated chemical contrast. 
(d,f) EELS spectra of points 0, 1, and 2, as marked in (c) and 
(e), respectively. 



substitutional Mn in the 3+ charge state. 

We determine the concentration ratio of the absorb- 
ing Mn 3+ ions for samples without and with Si by fitting 
a model constructed within the transfer matrix formal- 
ism of optical transmission and reflectivity, taking into 
account the doped (Ga,Mn)N layer, the undoped GaN 
buffer, and the sapphire substrate^ 

The contribution of the Mn 3+ ions to the dynamic di- 
electric function is modeled by damped Lorcntzian oscil- 
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FIG. 7: (Color online) Optical reflectivity spectra for 
(Ga,Mn)N without and with Si co-doping (upper and lower 
panel, respectively). The Mn concentration is 1.8% in both 
samples. The absorption feature at 1.41 eV is clearly resolved. 
The fitting results by the transfer matrix multilayer model are 
given by the dashed lines. 
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FIG. 8: (Color online) Optical transmission spectra for sam- 
ples with a Mn concentration x ~ 0.5%. Upper and lower 
graphs: data for samples without and with Si co-doping, re- 
spectively. Next to the main absorption line at 1.41 eV the 
phonon replicas of the line are observed. The fitting results 
by the multilayer model are given by the dashed lines. 



Label 


X 


E 


f N Mu3 + 


r 


Mn 2+ 




(%) 


(meV) 


(meV 2 ) 


(meV) 


/ ,-.90 / ^ \ 

(10 u / cm ) 


(Ga,Mn)N 


1.8 


1415 


13000 


6.1 




(Ga,Mn)N:Si 


1.8 


1414 


5300 


3.6 


4±1 


(Ga,Mn)N 


0.53 


1413 


950 


2.57 




(Ga,Mn)N:Si 


0.49 


1415 


490 


3.0 


0.6±0.1 



TABLE III: Values of parameters in Eq. 1 determined 
from reflectivity and transmission measurements for two se- 
ries of samples (the upper and lower panel, respectively). 
The concentration of Mn 2+ ions is calculated according to 

3 ,(Ga,Mn)N:Si_ a .(Ga,Mn)N r) where r ; g ^ ratk) Qf ^ ftf^^. 

values within a given series of samples. 



lators of the form. 
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Here / is proportional to the oscillator strength and T 
is the damping energy. The refractive index of GaN is 
modeled using the Scllmeier equation 3 ^ and the refrac- 
tive index of the sapphire is set at 1.8. The thicknesses 
of the GaN buffer and (Ga,Mn)N layer are adjusted to 
reproduce the observed Fabry-Perot interferences at a 
magnitude of thicknesses ratio as determined during the 
growth by in-situ ellipsometry. 

In Tabic Hill the fitted values of the parameters in Eq. 1 
are given. From the reduction of the Mn 3+ absorption 
we can evaluate the concentration of Mn 2+ ions assuming 
that electrons coming from Si shallow donors occupy the 
Mn 2+ /Mn 3+ midgap level. Under this assumption, the 
concentration of Mn 2+ ions reaches a level of 4x 10 20 per 
cm 3 for the highest employed flow of the Si precursor at 
the lowest Ga precursor flow. 

Results of ESR studies carried out in Linz and reported 
in Fig. M for (Ga,Mn)N and (Ga,Mn)N:Si with the Mn 
concentration x « 0.5% demonstrate the emergence of a 
characteristic Mn 2+ signal upon Si doping. In films with 
higher Mn concentrations, the line broadening, witness- 
ing the presence of Mn 2+ -Mn 3+ coupling, as discussed 
in Sec. IV, has hampered the detection of a Mn 2+ signal. 

A non-zero value of the orbital momentum, and the 
associated spin-orbit interaction specific to Mn 3+ ions, 
precludes their observation by ESR (Ref. |27|). At the 
same time, Mn 2+ ions, corresponding to orbital singlets, 
give rise to a specific ESR response! 27 ' 34 

The results of the optical and ESR studies are con- 
firmed also by XANES. In order to contribute to the 
determination of the Mn valence state, the position of 
the x-ray absorption edge and the pre-edge features have 
been considered as described in the XANES section of 
Ref. In particular here the issue of reduction of the 
charge state is addressed. As it can be appreciated in 
Fig. [lJJl a shift towards lower energies is visible for the Si 
co-doped sample, while the position of the pre-edge peaks 
remains unchanged. This demonstrates the fine calibra- 
tion of the incoming energy and suggests the presence of 
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FIG. 9: (Color online) Results of ESR studies at 2 K for 
(Ga,Mn)N without and with Si co-doping (lower and upper 
panel, respectively) showing the emergence of Mn 2+ signal on 
Si doping. The Mn concentrations is x fa 0.5%. 



some Mn in a valence state lower than 3+, possibly 2+. 
Due to the lack of model compounds that would allow 
to establish a precise relation between the edge shift and 
the valence state in nitrides, a quantitative statement 
cannot be given here. However, considering that in the 
case of 6-coordinated Mn ions (this example is taken as 
no data are reported for tetragonal Mn 3+ ) the ionic ra- 
dius of Mn 2+ is about 12% greater than the one of Mn 3+ 
(Ref. [35h and that no visible evolution of the Mn-N dis- 
tance is reported, we can expect a minority of Mn ions 
to be in the 2+ charge state. 
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FIG. 10: (Color online) Partial fluorescence yield XANES 
spectra (integrated K a i,2 fluorescence lines) for (Ga,Mn)N 
and (Ga,Mn)N:Si samples with the polarization vector paral- 
lel to the wurtzite c axis. Inset: focus on the shift of the main 
absorption edge between the two spectra. The position of the 
pre-edge peaks is constant. 

In conclusion, our SIMS, optical, XANES, and ESR 
studies show consistently that co-doping with Si increases 
the concentration of Mn ions in the 2+ charge state, 



which for the highest employed flow of the Si precur- 
sor (2 seem) and the lowest Ga precursor flow is up to 
4xl0 20 cm~ 3 , about 30% of the total Mn concentration 
for x = 3%. This evaluation substantiates the experimen- 
tal data presented in the next section. It is worth noting 
that a co-existence of Mn 3+ and Mn 2+ ions was detected 
also in x-ray magnetic circular dichroism in (Ga,Mn)N 
samples undoped with Sij2& pointing to the presence of 
residual or interfacial compensating donors. 



IV. MAGNETIC PROPERTIES 

According to our previous studies of (Ga,Mn)N with 
x < 1%,— the dependence of the magnetization M on 
temperature T, magnetic field H, and its orientation 
with respect to the wurtzite (wz) c-axis can accurately 
be described in terms of non-interacting Mn 3+ ions sub- 
stitutional of Ga. The good agreement between the ex- 
perimental data and the model confirms a weak com- 
pensation by residual impurities which, if present, would 
change the Mn charge state and thus the magnetic prop- 
erties. The Mn 3+ charge state is preserved in samples 
with higher Mn concentrations, where the persistence of 
a large anisotropy between the M(H) values at 1.85 K 
for the two sample orientations c _L H and c || H is 
evidenced in Figs. [TTTa) and ITTT b). However, a grad- 
ual enhancement of M(H)/M(b0 kOe) over the magni- 
tude expected for non-interacting spins is observed when 
increasing x up to 3% in both uniformly and digitally 
Mn-doped films, as seen in Figs. [TTTa). fTTTb). and [T2T a). 
These results demonstrate univocally that, in spite of the 
absence of band carriers, the dominant exchange interac- 
tion between Mn 3+ is ferromagnetic in (Ga,Mn)N. 

Interestingly, a rather different behavior is observed in 
the case of (Ga,Mn)N:Si, where the trapping of donor 
electrons changes the Mn charge from 3+ to 2+ and the 
spin state S from 2 to 5/2, for about 30% of the Mn ions 
at x rj 3%, as discussed in Sec. IIIC. According to the 
data collected in Figs. [TTtc) . fTTTd) . and[T27b) the increas- 
ing concentration of Mn 2+ ions results in the foreseen 
decrease of the magnetization anisotropy. Furthermore, 
as shown in Fig. IT2T b) , M (H) saturates slower than the- 
oretically anticipated for non-interacting Mn ions. This 
finding points to an antiferromagnetic character of the 
exchange coupling between Mn 2+ ions, and suggests that 
these ions may dominate in (Ga,Mn)N, when no ferro- 
magnetic interactions are detectedj 11 ' 37 

In the insets to Figs. ITTTa)-[rITd) the results of our 
search for the onset of a collective magnetic behavior 
in the samples with the highest Mn concentrations are 
given. A linear and ahysteretic M(H) dependence in 
weak magnetic fields is observed for both configurations, 
c±H and c\\H, pointing to the absence of spontaneous 
magnetization. These data imply that the ferromagnetic 
spin-spin couplings are too short-ranged to result in mag- 
netic ordering and, hence, in spontaneous magnetization 
at T > 1.85 K. 
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FIG. 11: (Color online) Magnetization at 1.85 K (normalized 
to its in-plane value at 50 kOe) as a function of the magnetic 
field applied perpendicular (circles) and parallel (triangles) 
to the wz c-axis of (Ga,Mn)N (a,b) and (Ga,Mn)N:Si (c,d) 
films where Mn is introduced either uniformly (a,c,d) or dig- 
itally (b). Solid and dashed lines are calculated according 
to the group theoretical model for non-interacting ions 
in wz-GaN for H perpendicular and parallel to the c-axis, 
respectively.— Insets: low-field magnetization loops. 



Quantitative information on J nn is gained here by ex- 
amining the dependence on the inverse temperature of 
the in-plane magnetic moment m{T) of Mn spins in 
GaN, as obtained by subtracting the value of m(T) mea- 
sured independently for a sapphire substrate (normal- 
ized by the corresponding sample weight). As reported 
in Fig. [H X(T) = M/H ~ 1/T for 150 < T ^ 350 K. 
This behavior indicates that the contribution to the mag- 
netic susceptibility from the (Ga,Mn)N films obeys the 
Curie law in this regime, x(T) = C/T. This depen- 
dence is expected if the spin pairs are either uncorrelatcd 
\Jnn\S 2 < fc B T, or strongly bound \ J nn \S 2 > k B T. 

In order to extract from these data the magnitude of 
J nn , we extend the previous model of a single substitu- 
tional Mn 3+ impurity in GaN (Ref. Gil) by considering 
a pair of Mn + ions coupled by an exchange interac- 
tion — JSiS2r^£2. the model discussed in details in the 
Appendix. Assuming a random distribution of Mn over 
the cation hep lattice, and allowing for the coupling be- 
tween nn spins we can evaluate M(T, H) at a given J nn 
and x. This approach implies, in particular, that for 
x = 3%, T < 350 K, and H = 1 kOe, X (T) = C/T if 
J nn > 10 meV. However, in this case, due to the propor- 
tionality of x(T) to the pair spin square, the magnitude of 
C is enhanced in comparison to the value Co correspond- 
ing to non-interacting spins. To evaluate experimentally 
Cnorm = C/Co, we consider that its magnitude can be 
determined from the magnetic moment m(T, 1 kOe) mea- 
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FIG. 12: Comparison of magnetization saturation as a func- 
tion of the magnetic field between (a): (Ga,Mn)N - no anti- 
ferromagnetic interactions and (b): (Ga,Mn)N:Si - with an- 
tifcrromagnetic interactions. The relative experimental error 
is about one half of the point size. Solid lines and the dashed 
line in (b) are calculated for non-interacting Mn 3+ and Mn 2+ 
ions, respectively. 
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FIG. 13: (Color online) Magnetic susceptibility M/H for 
c _L H as a function of the inverse temperature for uniformly 
(solid symbols) and digitally (open symbols) doped (Ga,Mn)N 
(circles) and (Ga,Mn)N:Si (triangles) as well as GaN/sapphire 
control sample (light grey circles). Solid lines serve to deter- 
mine the Curie constants at high temperatures. Deviations 
of their values form zero at 1/T = show the accuracy of the 
substrate subtraction. 



sured in-plane without knowing the exact value of the 
volume occupied by the Mn spins, if the magnitude of 
the in-plane m(1.85 K, 50kOe) is employed to obtain the 
Mn content x — and thus Co — for particular samples. 
Following the outcome of the experimental results for 
(Ga,Mn)N:Si (Sec. IIIC), demonstrating the presence of 
Mn 2+ ions, their relative contribution to M(T, H) is de- 
termined from the magnitude of the magnetic anisotropy. 
As summarized in Fig. 1141 C norm > 1 for all studied 
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FIG. 14: (Color online) Normalized Curie constant C n0 rm as 
a function of the Mn content for uniformly (solid symbols) 
and digitally (open symbols) doped (Ga,Mn)N (circles) and 
(Ga,Mn)N:Si (triangles). The solid lines are computed as- 
suming a random distribution of Mn 3+ ions, and either fer- 
romagnetic or antiferromagnetic strong coupling between the 
nearest neighbor (nn) Mn spins, | J nn |S 2 3> ksT. Dashed and 
dotted lines are calculated assuming that 4 • 10 20 cm -3 Mn 
ions are in the 2+ charge state, and nn interactions are fer- 
romagnetic except for nn Mn 2+ pairs for which J nn = or 
— JnnS 2 ^> ksT, respectively. 

samples, hinting to the presence of a considerable fer- 
romagnetic spin-spin interaction. The theory presented 
in the Appendix describes quite well the magnitude of 
C noxm (x) for the uniformly doped (Ga,Mn)N films, point- 
ing to Jnn > 10 meV, in general agreement with the 
results of the ab initio studies outlined in Fig. Q] Fur- 
thermore, a low-temperature upturn of the experimental 
points over the C/T dependence, visible in Fig. [TBI below 
~100 K, suggests the existence of an additional weak 
ferromagnetic coupling between more distant neighbors. 
Moreover, the experimental values of C norm (x) in the case 
of digital <5Mn-doping are higher than theoretically ex- 
pected. 

In order to clarify the different magnitude of C n0 rm in 
uniformly and digitally doped films, we refer to Sec. IIIB, 
where the detailed investigation of the Mn distribution 
for the two samples with the highest Mn concentra- 
tion, respectively (Ga,Mn)N (x = 3.1%) and the digi- 
tal (Ga,JMn)N (x av = 2.6%), have been shown. The 
data demonstrate the presence of a spatially modulated 
Mn concentration in the digitally Mn-doped films. Due 
to a non-linear dependence of the Curie constant on the 
Mn concentration in the presence of ferromagnetic in- 
teractions, such a non-random distribution of Mn ions 
increases the apparent value of C n0 rm, particularly if the 
system is close to a ferromagnetic instability. This in- 
terpretation is supported by a much smaller effect in the 
films with lower values of cc av , and thus far from the fer- 
romagnetic instability. 

Finally, we comment on the magnitudes of C noTm (x) 
in Si doped samples. Here, we have ferromagnctically 



coupled Mn 3+ -Mn 3+ and Mn 3+ -Mn 2+ pairs as well as 
antiferromagnetically interacting Mn 2+ -Mn 2+ pairs. As 
shown in Fig. [TH the theory developed for such a case 
and summarized in the Appendix is consistent with the 
data for Si-doped samples. 



V. CONCLUSIONS 

In this work, we have verified experimentally the pres- 
ence of a strong ferromagnetic coupling between neigh- 
boring Mn spins in (Ga,Mn)N, supporting a very signif- 
icant number of ab initio studies ££iL Since the Mott- 
Hubbard localization precludes carrier hopping between 
magnetic ions, a ferromagnetic super-exchange consti- 
tutes the relevant microscopic coupling mechanism^ 
However, according to our findings, the range of this in- 
teraction is too short to produce a long-range ferromag- 
netic ordering - at least above 1.85 K - in samples with 
3% of randomly distributed substitutional Mn cations. 
Co-doping with Si may a priori result in a ferromag- 
netic double exchange, but apparently Andcrson-Mott 
localization in the Mn impurity band renders this mech- 
anism inefficient in the range of Mn contents explored 
so far by us. If, owing to a large density of donor-like 
defects or impurities the concentration of Mn 2+ prevails, 
antiferromagnetic super-exchange becomes the dominant 
spin-spin coupling mechanism. This situation has pre- 
sumably taken place in recently studied (Ga,Mn)N films 
with x up to 36% (Ref. [Til ) and a time ago in the case of 
Ini-^Mn^As layers with x up to 18%4I 
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Appendix: Theoretical evaluation of the Curie 
constant 

We evaluate the Curie constant C norm for a random 
distribution of the Mn 3+ and Mn 2+ ions. While our 
model can be applied for a general situation, we dis- 
cuss the case describing our data in the relevant temper- 
ature range 150 < T ^ 350 K, i.e., we assume that the 
nearest neighbour (nn) Mn 3+ -Mn 3+ and Mn 3+ -Mn 2+ 
pairs form ferromagnctically oriented dimcrs, whereas nn 
Mn 2+ -Mn 2+ pairs are uncoupled. 

In the paramagnetic region well above the ordering 
temperature, (high temperature limit) the magnetic sus- 
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ccptibility is expected to obey the Curie- Weiss law, 

C 



X 



T-6 r 



Cn = N 



( gf i B ys(s + i) 

3k B 



(5.2) 
(5.3) 



where C and 6c are the Curie constant and Curie- Weiss 
temperature, and N is the concentration of magnetic 
ions with spin S. However, in random magnetic alloys, 
where the interactions between spin pairs show a large 
dispersion owing to strong variations of the spin-spin dis- 
tances, the magnitudes of C and 9c may depend on the 
temperature i 38 ! 39 ' 42 

We consider the case of dilute magnetic semiconduc- 
tors (DMSs) and dilute magnetic oxides (DMOs). In the 
absence of band carriers that could mediate long range 
spin-spin interactions, the strength of the exchange cou- 
plings decays rather fast with the spin-spin distance. In 
such a case, the exchange between magnetic ions occupy- 
ing the nn cation positions dominates. Accordingly, be- 
low the percolation limit for the nn interaction (x < 18% 
for fee and hep lattices), the magnetic response can be 
evaluated as a sum of contributions coming from various 
types of clusters: isolated spins, nn pairs, nn triads, 
whose relative importance depends on The magnetic 
response of such small clusters can be easily calculated 
for given values of the exchange integral J, m , tempera- 
ture T, and magnetic field H. Possible interactions be- 
tween more distant pairs of magnetic ions are neglected 
within this approach. Qualitatively, at low temperatures, 
k^T -c |</nn|j at which the nn pairs arc locked, the Curie 
constant is reduced (C < Co) or enhanced (C > Co) for 
antiferromagnetic interactions (J nn < 0) and ferromag- 
netic interactions (J nn > 0), respectively. 

We are interested in the case of (Ga,Mn)N and 
(Ga,Mn)N:Si where, in general, both Mn 2+ (S = 5/2) 
and Mn 3+ (S = 2) ions are present with the concen- 
trations x Mn 2+ and a; Mn 3+, respectively (x = x Mn 2+ + 
x Mn s+). Here, we have ferromagnetically coupled Mn 3+ - 
Mn 3+ and Mn 3+ -Mn 2+ pairs as well as antiferromag- 
netically interacting Mn 2+ -Mn 2+ pairs. To describe this 
case, we generalize the approach put previously forward 
for II- VI DMSs^ The probability that a given spin, e.g., 
a Mn 3+ ion, is in a cluster belonging to the configuration 
r cl is given by, 

Pr cl = n rcl (x Mn2+ )"M„2 + (x Mn3+ )"M„3 + -i (1 _ x)VTc1 ^ 

(5.4) 



where n Mn 2+ and n Mn 3+ are the numbers of manganese 
ions in 2+ and 3+ charge states belonging to this clus- 
ter and n rcl , v Tcl are parameters taken from Refs. l39l and 
|43| . Then one should take into account all possible com- 
binations of arrangements of Mn 2+ ions within r c i and 
calculate the total spin S c \ corresponding to the ground 
state (assuming that |J nn 3> k^T). In the case of van- 
ishing interactions between Mn 2+ ions we have to com- 
pute also the effective cluster size n c \ (some of the ions 
may be disconnected from the initial cluster because of 
J Mn 2+_ Mn 2+ = ). In this way we obtain the probabil- 
ity matrix Ps i,n i from which the Curie constant can be 
calculated, 



C = xN 



3fc B 



E 

S c i,n c 



Ps al ,n c iS c l(S c \ + 1) 



rid 



(5.5) 



where Nq is the cation concentration. 

In order to obtain the magnetization M (T, H) of 
(Ga,Mn)N in the presence of interacting nn magnetic 
centers, we extend the previous model of a single sub- 
stitutional Mn 3 + impurity in GaN (Refs. [T3IlMi3) by 
considering a pair of Mn 3+ ions coupled by an exchange 
interaction H (12) = —J nn SiS2- Then, the energy struc- 
ture of such a pair can be described by the Hamiltonian 



H = H(1) + H(2) + H(12), 



(5.6) 



where H(i) (i = 1,2) accounts for the single Mn 3+ 
ion (L = 2, S = 2) in GaN with the trigonal crystal 
field of the wurtzite structure and the Jahn- Teller dis- 
tortion taken into account (for details, see Ref. [13). In 
the coupling scheme employed, the base states for the 
pair are characterized by the set of quantum numbers 
, ms 1 , TiL 2 , Tis 2 ) ■ The energy level scheme is calcu- 
lated by a numerical diagonalization of the full 625 x 625 
Hamiltonian (|5.6[) matrix, which allows to obtain an av- 
erage magnetic moment of the Mn ion belonging to the 
pair (m pa j rs ). Assuming a random distribution of Mn 
ions over the cation sites (the hep lattice) and allowing 
for the coupling between the nn spins we can then eval- 
uate M(T, H) at given J nn and x, 



M = HBxN [(m sing i es )P ncl=1 + (mpairs) (1 - Pn ol =l) 



(5.7) 
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